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Abstract The anomalous X-ray pulsar 4U 0142+61 was 
recently detected in the mid infrared bands with the 
SPITZER Observatory (Wang, Chakrabarty & Kaplan 
2006). This observation is the first instance for a disk 
around an AXP. From a reanalysis of optical and in- 
frared data, we show that the observations indicate that 
the disk is likely to be an active disk rather than a pas- 
sive dust disk beyond the light cylinder, as proposed in 
the discovering paper. Furthermore, we show that the 
irradiated accretion disk model can also account for all 
the optical and infrared observations of the anomalous 
X-ray pulsars in the persistent state. 
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1 Introduction 

Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray 
Repeaters (SGRs) (Mereghetti et al. 2002; Hurley 2000; 
Woods & Thompson 2006) are characterized mainly by 
their X-ray luminosities (L x ~ 10 34 — 10 36 erg s _1 ), 
which are orders of magnitude higher than their rota- 
tional powers E ro t = Ififi. Their spin periods are clus- 
tered to a very narrow range (5 — 12 s). Most of SGRs 
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and AXPs show repetitive, short (^ls) super-Eddington 
bursts with luminosities up to 10 erg s _1 . Three giant 
flares with peak luminosities L p > 10 44 erg s _1 and dura- 
tions of a few minutes were observed from three different 
SGRs (Mazets et al. 1979; Mazets et al. 1999; Hurley et 
al. 1999; Palmer et al. 2005). 

Energetics of the bursts strongly indicate a magnetar 
mechanism for these bursts. Magnetar models (Thomp- 
son & Duncan 1995; Thompson & Duncan 1996) adopt 
strong magnetic fields with magnitude _B» > 10 14 G on 
the stellar surface to explain the burst energetics. These 
bursts are likely to be local events near the neutron star 
and their energies could be originating either from the 
dipole component or from the higher multipoles of the 
stellar magnetic field. In the current magnetar models, 
it is the dipole component of the magnetic field which 
must be of magnetar strength to account for the spin 
down properties of the AXPs and SGRs. 

In the alternative fallback disk model (Chatterjee, 
Hernquist, & Narayan 2000; Alpar 2001), the source of 
the X-rays is accretion onto the neutron star, while the 
optical/IR light originates from the accretion disk. Rota- 
tional evolution of the neutron star is determined by the 
interaction between the disk and the dipole component 
of the neutron star's magnetic field (B* ~ 10 12 — 10 13 
G). Fallback disk models can account for the period 
clustering of AXPs and SGRs as the natural outcome 
of disk-magnetosphere interaction during their lifetimes 
(Alpar 2001; Ek§i & Alpar 2003). As higher multipole 
fields rapidly decrease with increasing radial distance (as 
r~ 5 for quadrupole component), the dipole component 
of the magnetic field determines the interaction and the 
angular momentum transfer between the disk and the 
neutron star. The strength of magnetic dipole field in 
order to explain the period clustering of the AXPs and 
SGRs over their M history is B* ~ 10 12 - 10 13 G (Alpar 
2001; Ek§i & Alpar 2003). Therefore, disk models are 
consistent with the magnetar fields in the quadrupole or 
higher components, while they are not compatible with 
magnetar fields in the dipole component. 
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Fallback disk models can also explain the enhance- 
ments observed in the persistent luminosities of SGRs 
and AXPs. The X-ray enhancement of the SGR 1900+14 
following its giant flare can be explained by the relax- 
ation of a disk which has been pushed back by a pre- 
ceding burst (Ertan & Alpar 2003). The same model 
with similar disk parameters can also reproduce the cor- 
related X-ray and IR enhancement of AXP 2259+58, 
which lasted for ~ 1.5 years, if this is triggered by a 
burst, with a burst energy estimated to have remained 
under the detection limits (Ertan, G6gu§ & Alpar 2006). 
The suggestion of fallback disks has motivated observa- 
tional searches for disk emission in the optical and IR 
bands, and resulted in various constraints on the mod- 
els. Some of the AXPs were observed in more than one 
IR band (Hullcman ct al. 2001; Israel ct al. 2002; Wang 
& Chakrabarty 2002; Kaspi et al. 2003; Israel et al. 2003; 
Hulleman et al. 2004; Israel et al. 2004; Tarn et al. 2004; 
Morii ct al. 2005; Durant & van Kerkwijk 2006a). AXP 
4U 0142+61 is the source with most extended obser- 
vations, as it was also observed in the optical R and 
V bands (Hulleman et al. 2000; Hulleman et al. 2004; 
Dhillon et al. 2005, Wang et al. 2006). Discovery of mod- 
ulation in the R band luminosity of 4U 0142+61 at the 
neutron star's rotation period P=8.7 s, with a pulsed 
fraction 27 % (Kern & Mertin 2002; Dhillon et al. 2005), 
is particularly significant. This fraction is much higher 
than the pulsed fraction of the X-ray luminosity of this 
source, indicating that the origin of the pulsed optical 
emission cannot be the reprocessed X-rays by the disk. 
Magnetospheric models for these pulsations can be built 
with either a dipole magnetar field or within disk-star 
dynamo model (Cheng & Ruderman 1991), in which a 
magnetospheric pulsar activity is sustained by a stellar 
dipole field of ~ 10 12 G and a disk protruding within 
the magnetosphere. Ertan & Cheng (2004) showed that 
this pulsed optical component of the AXP 4U 0142+61 
can be explained by both types of magnetospheric mod- 
els. Thus the presence of strong optical pulsations from 
the magnetosphere does not rule out the possibility of a 
fallback disk together with 10 12 — 10 13 G surface dipole 
magnetic field. 

In the present work, we concentrate on the unpulsed 
optical/IR emission from the AXPs and SGRs in their 
persistent states, and test the expectations of the irra- 
diated accretion disk model through the observations in 
different optical/IR energy bands (V, R, I, J, H, K, and 
K s ). The optical/IR emission expected from the irra- 
diated fallback disks was first computed and discussed 
by Pcrna, Hcrnquist & Narayan (2000) and Hulleman, 
van Kerkwijk & Kulkarni (2000). Using similar irradia- 
tion strengths, Perna et al. (2000) and Hullcman et al. 
(2000) found similar optical fluxes that remain well be- 
yond those indicated by the observations of AXP 4U 
0142+61 and AXP IE 2259+586. To explain this re- 
sult, Perna et al. (2000) suggested that the inner disk 
regions could be cut by an advection dominated flow, 



while Hullcman et al. (2000) concluded that the then 
existing optical data of the AXP 4U 0142+61 (in I,R,V 
bands) can only be accounted for by an extremely small 
outer disk radius, around a few x 10 9 cm (see Ertan et al. 
2007 for detailed discussion). In present work, we show 
that the optical/IR data of the AXPs can be explained by 
the irradiated accretion disk model without any implau- 
sible constraints on the outer and inner disk radii. The 
main reason for the difference between our results and 
those of earlier works is that both Hulleman ct al. (2000) 
and Perna et al. (2000) assumed a particular irradiation 
strength, while we keep it as a free parameter considering 
the observations of the low mass X-ray binaries (LMXBs) 
which indicate varying irradiation strengths (see Ertan 
& Cah§kan (2006) for a detailed discussion) We give the 
details of the disk model in § 2. We discuss our results 
in § 3, and summarize the conclusions in § 4. 



2 Optical/IR Emission from the Irradiated Disk 

Model fits to the contemporaneous X-ray and IR data 
of the AXP IE 2259+586 favor the irradiated disk mod- 
els, though they do not exclude the nonirradiated thin 
disk model (Ertan, Gogiis, & Alpar 2006). We start by 
assuming that the AXP disks are irradiated and include 
the irradiation strength as a free parameter through our 
calculations. 

We calculate the disk blackbody emission taking both 
the intrinsic dissipation and the irradiation flux into ac- 
count. A steady disk model is a good approximation for 
the present evolution of the AXP and SGR disks in their 
persistent states. For a steady thin disk, the intrinsic 
dissipation can be written as 
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(see e.g. Frank ct al. 2002) where M is the disk mass 
flow rate, M is the mass of the neutron star and R is the 
radial distance from the neutron star. In the absence of 
irradiation, the effective temperature T ff of the disk is 
proportional to i?~ 3 / 4 for a given M. For an irradiated 
disk, the irradiation flux can be written as 
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(Shakura & Sunyaev 1973), where r\ is the conversion 
efficiency of the rest mass energy into X-rays, e is the 
X-ray albedo of the disk face, c is the speed of light and 
H is the pressure scale height of the disk. Since the ratio 
H/R is roughly constant along the disk, the irradiation 
strength is expected to be constant along the disk. De- 
pending on the geometry and the temperature of the 
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scattering source, efficiency of the irradiation could vary 
from source to source. Estimates for the parameter C are 
usually in the range 1CP 4 — 10~ 3 for the low mass X-ray 
binaries ( de Jong at al. 1996; Dubus et al. 1999; Ertan 
& Alpar 2002). 

For a constant C along the disk, the irradiation tem- 
perature Tj rr = (Fi^/a) 1 / 4 is proportional to i? -1 / 2 (Eq. 
2). For small radii, dissipation is the dominant source 
of the disk emission. At a critical radius R Cl the irra- 
diation flux becomes equal to the dissipation rate, and 
beyond R c , the dominant heating mechanism is X-ray ir- 
radiation. Equating F irr to D (Eqs. 1 and 2), the critical 
radius is found to be 



3 GM, 
2~C^ 



10" 



C 



3 x 10 9 cm. 



(4) 



The effective temperature profile of the disk can be writ- 
ten as 

T«-(° + *±)"' (5) 
a a 



where a is the Stefan-Boltzmann constant. 

To find the model disk flux in a given observational 
band, we integrate the calculated blackbody emissions 
of all radial grids radiating in this band. For comparison 
with data, we calculate the model disk fluxes along the 
optical/IR bands V, R, H, I, J, K and K s . For all sources 
we set cos i = 1 where i is the angle between the disk 
normal and the line of sight of the observer. We equate 
the disk mass flow rate M to the accretion rate onto 
the neutron star, thus assuming the mass loss due to 
the propeller effect is negligible. We first adjust M to 
obtain the observed X-ray flux. Next, using this value 
of M and taking the strength of the magnetic dipolc 
field B* = 10 12 G on the surface of the neutron star we 
calculate the Alfven radius Ra which we take to be the 
inner radius of the disk. Then, we look for a good fit 
to the overall available optical/IR data by adjusting the 
irradiation strength C within the uncertainties discussed 
in § 3. For the sources with more than one detections 
in a particular IR/optical band we take the observation 
nearest to the date of X-ray observation. 



3 Results and Discussion 

We summarize our results in Table 1. For each source, 
the first column gives the unabsorbed flux data obtained 
from the observed magnitudes and the estimated Ay val- 
ues (for references see Table 1), and the second column 
gives the model fluxes. For the AXP 4U 0142+61, rea- 
sonable range of reddening is 2.6 < Ay < 5.1 (Hulleman 
et al. 2004). We obtain a good fit with A v = 3.5 (Fig. 
1). Table 1 shows that the irradiated steady disk model 
is also in agreement with all the other AXPs observed in 
the optical and IR bands. The parameters of the model 
for each source are given in Table 2. 



At present, AXP 4U 0142+61, which has been ob- 
served in nine different optical/IR bands from 8 ^tm to 
V in the same X-ray luminosity regime, seems to be the 
best source to study the properties of AXPs in the per- 
sistent state. The irradiation parameter C obtained from 
our model fits remains in the range (10~ 4 < C < 10~ 3 ) 
estimated from the observations of LMXBs and the disk 
stability analyses of the soft X-ray transients (de Jong et 
al. 1996; Dubus et al. 1999; Ertan & Alpar 2002). Within 
the critical radius R c given by Eq. 4, dissipation is the 
dominant heating mechanism. For the disk model of the 
AXP 4U 0142+61, R c ~ 3 x 10 9 cm and R in = 1 x 10 9 
cm. The innermost disk emitting mostly in the UV bands 
also contributes to the optical emission. The radial dis- 
tance at which the disk blackbody temperatures peak 
at the optical bands (R,V) is about 10 10 cm. Peak tem- 
peratures of the IR bands from I to SPITZER mid IR 
bands (4.5 fxm and 8 /im) lie between fi~2x 10 10 cm 
and R <~ 10 12 cm. There are several observations of this 
source in some of the IR bands; we adopt the data near- 
est to the X-ray observation epoch used in our analysis 
here. 

For AXP J1708-40, Durant and van Kerkwijk (2006a) 
recently found that the previously reported IR data in 
K s , H, and J bands are likely to be a background star. 
They found another object within the positional error 
cycle and argued that this second object is more likely 
to be the IR counterpart to the AXP J1708-40. For this 
source, we adopt the IR (K s , H, J) data set reported by 
Durant and van Kerkwijk (2006a). 

For the the AXP IE 2259+586, we use the data taken 
before the X-ray enhancement phase of this source (Hulle- 
man et al. 2001). This source was detected in K s band 
and there are upper limits for I and R bands. Our model 
flux values are three and ten times below the upper limits 
reported for I and R bands respectively. 

AXP IE 1048-59 was detected in K s , H and I bands 
(Wang & Chakrabarty 2002). Observed X-ray flux from 
this source between December 2000 to January 2003 show 
a variation within a factor of 5 (Mereghetti et al. 2004). 
We use the X-ray flux obtained from the nearest X-ray 
observation to the date of the IR observations when the 
source was in the persistent state. 

AXP IE 1841 was detected only in the K s band, and 
there is a high upper limit in the R band (Wachter et 
al. 2004). Model estimates in other optical/IR bands for 
this source (and the other AXPs) can be tested by future 
optical and IR observations. 

We note that our model does not address the uncer- 
tainties in the innermost disk interacting with the mag- 
netosphcrc, the contributions from the magnetosphcric 
pulsed emission which is known to have a fraction about 
27 % in the R band, and the shielding effects which might 
decrease the irradiation strength at some regions of the 
disk. All of our optical and IR flux values remain within 
30 % of the data points, which is a reasonable fit consid- 
ering the uncertainties mentioned above. 
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Fig. 1 Energy flux data and irradiated disk model values for the AXP 4U 0142+61 in the optical and infrared bands (B, V, 
R, I, J, H, Kg, 4.5 /ttm and 8 fim). Filled circles are the observed (absorbed) data (taken from Hulleman et al. 2000 (V, R, 
I), Hulleman et al. 2004 (B,Kg), Morii et al. 2005 (H, J), Wang et al. 2006 (4.5 /im and 8 /mi), and triangles are dereddened 
data with Ay — 3.5. Squares are the irradiated disk model energy flux values (see § 2 for details). 



Table 1 The data flux values were calculated by using the magnitudes and ^4v values given in the references below. 
REFERENCES: (J1708-40) Durant & van Kerkwijk 2006a, Rea et al. 2003; (IE 2259+586) Hulleman et al. 2001, Woods et 
al. 2004; (IE 1841-45) Wachter et al. 2004, Morii et al. 2003; (IE 1048-59) Wang & Chakrabarty 2002, Mereghetti et al. 2004 
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Table 2 The parameters of the irradiated disk model which gives the optical/IR flux values seen in Table 1. For all the 
sources, we set cos i = 1 where i is the inclination angle between the disk normal and the line of sight of the observer, and 
we take the outer disk radius R ou t = 2 x 10 12 cm. See Section 3 for details 
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Since there is no detections in short wavelength opti- 
cal bands for the AXPs (except for AXP 4U 0142+61), 
model fits are not sensitive to the chosen inner disk radii. 
For the model fits, we equate the inner disk radii to the 
Alfven radii (Table 2) corresponding to a magnetic field 
with magnitude £?* = 10 12 G on the stellar surface and 
the accretion rates derived from the estimated X-ray 
luminosities (see Table 1 for references). On the other 
hand, optical data of the AXP 4U 0142+61 in R and V 
bands provide a constraint for the inner disk radius, and 
thereby for the strength of the magnetic dipole field of 
this source (see § 4). 



4 Conclusion 

We have shown that the optical, infrared and X-ray ob- 
servations of the AXPs in their persistent states can be 
explained with irradiated disk models. Among the AXPs, 
4U 0142+61 is currently the only source which provides 
an upper limit for the inner disk radius through its op- 
tical (R,V) data. For the best model fit for this source, 
which we have obtained with Ay = 3.5, the model in- 
ner disk radius (~ 10 9 cm) is around the Alfven radius 
for the accretion rate, estimated from the X-ray lumi- 
nosity, together with a dipole magnetic field strength 

~ 10 12 G on the neutron star surface. Nevertheless, 
it is possible to obtain reasonable fits by increasing the 
inner disk radius and decreasing the reddening accord- 
ingly. For Ay = 2.6, the minimum value of the reddening 
in the reasonable range 2.6 < Ay < 5.1 (Hulleman et al. 
2004), we obtain the best fit with R m ~ 8 x 10 9 cm 
which corresponds to the maximum reasonable dipole 
field strength S, ~ 4 x 10 13 G on the pole and half 
of this on the equator of the neutron star. Nevertheless, 
we note that these limits could be increased depending 
on the amount of possible mass loss due to propeller 
effect and/or on how much the inner disk radius pen- 
etrates inside the Alfven radius (see Ertan et al. 2007 
for a detailed discussion for 4U 0142+61). On the other 
hand, even including these possibilities, very recent anal- 
ysis concluding Ay = 3.5 ± 0.4 (Durant & van Kerkwijk 
2006b) implies surface dipole magnetic field strengths 
less than about 10 13 G. The magnetar fields (£?* > 10 14 
G) in multipoles, that could be responsible for the burst 
energies, are compatible with this picture, while the op- 
tical data excludes a hybrid model involving a disk sur- 
rounding a dipole field of magnetar strength. 

On the other hand, existing IR data of the AXPs, in- 
cluding recent observations of 4U 0142+61 by SPITZER 
in 4.5 ^m and 8 fim bands (Wang et al. 2006), do not put 
an upper limit for the extension of the outer disk radius 
R out - The lower limit for R out provided by the longest 
wavelength IR data of the AXP 4U 0142+61 is around 
10 12 G. Further observations in the longer wavelength 
infrared bands by SPITZER space telescope will provide 



valuable information about the structure and possibly 
the extension of the fallback disks around these systems. 
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